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ABSTRACT 


This thesis reports on a study involving the simula- 
tion and control of a shipboard generating system 
(l10O00KW). The first part of the paper deals with the 
Simulation of the existing system which includes a 
steam turbine, synchronous generator, voltage regulator 
and governor. Data on the system was obtained from the 
various manufacturers involved. Several simulations 
were made using different ways of modeling a synchronous 
machine that have been proposed in the literature, and 
their computational requirements compared. The models 
were comparable in the level of physical detail being 
moeeosceneeod ana al’ yge~ro fernd to acres crise wlosel 
with test results. 

The second part of the thesis reports on an attempt 
to control the system using a digital computer directly 
in the control loops. The simulation was repeated with 
the turbine and generator modeled as indicated above, 
while the voltage regulator and governor functions 
were implemented using relatively simple algorithms. 
For the degree of control needed in the shipboard 
environment these algorithms proved to be reasonably 


efficient and to give good results. 
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fe LN LRODUCTION 


The modern steam powered ship is totally depend- 
ant on its ability to generate electric power. Steam 
powered may even be a misnomer since without electric 
power steam cannot be generated. Presently all that 
is steam powered is the turbogenerator, the main feed 
pumps, and the heating system. In the engineering plant 
the fuel oil pumps, the feed booster pumps, the condens- 
ate pumps, the fire and flushing pumps, and the cooling 
water pumps are electric powered. Even the automatic 
control system for the boiler requires electric power 
and cannot use steam power directly. Outside the 
engineering spaces everything is powered by electricity 
from the large electric motors for the steering system, 
for the gun and missile loaders and launchers, and for 
the anchor windless; to the complicated electronic 
communication, search and surveillance devices. With- 
out a reliable source of electric power a modern ship 
is useless. 

In order to reduce the operating costs of ships, 
automation is expanding into the shipboard environ- 
ment. One area being investigated is the use of the 
digital computer to reduce the number of personnel 
required to operate the ship's engineering plant. The 
possibilities range from just using the computer for 


the collection and logging of data all the way to the 





computer controlling the entire engineering plant. If 
such a computer were added to the engineering plant 

how could it best improve the generation of electricity 
aboard ship? One possibility is to replace the analogue 
controllers of the turbogenerators with a computer 
controller. This appears promising because of the special 
properties of a shipboard electrical system. On board 
ship the turbogenerators are located adjacent to each 
other. Normally the switchboard is located in an adjacent 
compartment with each generator connected to the bus line 
in this switchboard. The various shipboard loads are 
powered from this bus line. The generators are either 
totally connected at this point or the bus line is broken 
ana each generatuLr Carries a seperate part of tne i1oad 
with no interconnection. In either case, mathematically 
the electric system would appear to be a single machine 
connected to a single load. In the shipboard system the 
Main engines use about ninty-five percent of the steam 
produced. This means that the turbogenerators have such 

a large supply of steam that for analysis purposes it 
appears to be infinite. Thus, the boiler response does 
not enter into the turbogenerator system. Finally and 
most important, onboard ship the ratio of the change in 
electrical power required to the available electrical 
power is many times larger than is normally found in a 


large utility system. This results in larger fluctuations 





in the turbogenerator parameters of frequency and terminal 
voltage. Thus the action of the controllers of these 
parameters are more easily seen and the extreme accuracy 
required ina large utility system is not expected to be 
met by the shipboard controllers. 

In a shipboard environment there should be a simple 
efficient digital control scheme for the turbogenerator 
parameters of frequency and terminal voltage which would 
meet or exceed the response of the present controllers. 
Currently the frequency is controlled by a Watt flyball 
governor detecting speed changes and a hydraulic servo 
system translating these speed changes into changes in 
throttle position. This changes the power of the turbine 
Sc. FSCUlNS Cae Spec we Sa.ed Vaiuc. THE cCuDeent speed 
controller maintains speed within about four per cent of 
rated value under normal load conditions. It returns 
the speed to within two per cent of rated value in less 
than two seconds. The terminal voltage is controlled 
by a zener diode voltage sensing bridge, and a current 
transformer. The generator has a static excitor which 
uses the output of the current transformer as the input 
to a magnetic amplifier. The voltage sensing bridge 
controls the gain of the magnetic amplifier. The current 
transformer provides a vector sum of the terminal voltage 
and current, which aids the control fUnCeE Lon. Together 
the voltage controller maintains voltage within one per 


Senc uncer normal conditions. 





To find and test such a digital control scheme 
a model of the turbogenerator will be developed. A 
typical shipboard turbogenerator will be used for the 
basis of the model. The complete model will be tested 
against actual test results. A digital control scheme 
will be proposed. The digital controllers will replace 
the current controller in the model, and more simula- 
tions will be tested. The responses of the controllers 
will be compared. 

During the development of the generator model 
the efficiency of the use of Parks' Transformation was 
questioned. To determine the order of magnitude of the 
Savings and where tnese savings came from it was decided 
to try to simulate the mndel both with and withnnt +he 


use of Parks' Transformation. 





mes toe DEVEZOPMENT OF THE MODEL 
2.-I BACKGROUND 

To insure authenticity of the control scheme 
the model is based on an actual shipboard type turbo- 
generator. The test results used to determine machine 
constants were provided by Allis Chalmers from the 
acceptance tests for the turbogenerator installed 
at Great Lakes Naval Station Training Plant. This 
is a mock up of a typical shipboard engineroom and is 
totally operational. The same model turbogenerator may 
be found in the Leahy Class Guided Missile Teigewes 
(DLG-16) . 

The turbogenerator consists of a DeLaval 1442 BHP 
horizontal shaft condensina turbine. DeLaval speed 
reducing gears, an Allis Chalmers 1000 KW 450 Volt, 60 
Cycle generator and static excitor, and Allis Chalmers 
auxiliary equipment. 

Figure 1 shows the major parts of the shipboard 
environment. There are four boilers, two main engines, 
three turbogenerators, and much more major and minor 
auxiliary equipment, [1]. The boilers produce 
approximately 120,000 pounds of superheated steam (1200 
PSI and 940°F) per hour. The majority of the steam goes 
directly to the main engines. A very small part is used 
for the few pieces of steam powered auxiliary equipment, 


primarily the main feed pumps, and a much smaller 
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Guantity 1s reduced to 150 PSI for hotel services such 

as heating, cooking, and hot water. Of course there is 
approximately 9000 pounds of steam per hour that is used 
to generate electricity. The turbogenerators receive 
steam via a reducer, thus, being supplied at 440 PSI and 
740°%F. The electrical output of the generators goes 
directly to the main switchboard where all the generators 
are connected to the one bus line. 

From this information two basic assumptions will be 
drawn. First since the turbogenerators use such a small 
portion of the available steam and the steam that they use 
1s provided via a reducer, the turbogenerators appear to 
have an infinite supply of steam at constant temperature 
and pressure. Second since the turbogenerators are 
physically and electrically close they appear as a single 
machine supplying a single load. 

Throughout the paper the quantities will be expressed 
in per unit notation unless otherwise specified. The bases 


for the per unit system are as follows: 


Volt Amperes IT Oak = 1250 KVA 

Volts 1 p.u. = 260 Volts 

Amperes 1 p.u. = 1604 Amps 

Rotor Volts 1 p.u. = 54065 Volts 

Rotor Amperes sl plu. = £1.66 Amps 

Frequency 1 p.u. = 377 Radians per second 
Steam Flow 1 p.u. = 10,000 lbs per hour 
Impedence ie... = 2.1621 Ohms 





Zac ee, TURBINE 

The DeLaval turbine of the turbogenerator is a 9 
stage 10,009 RPM condensing turbine [2]. It receives 
steam at 440 PSI and 740°F and exhausts to 4 inches of 
mercury absolute. (Condensors operate at pressures below 
atmospheric pressure and thus, to more accurately define 
the vapor point in the condensor, the pressure is measured 
in inches of mercury.) The turbine is capable of carry- 
ing its rated load continuously. The system is shown in figure 2. 

Steam flow to the turbine is controlled by the 
throttle admitting steam to three nozzles progressively. 
In periods of low demand the steam is admitted only to 
the first nozzle. As the power required increases 
steam is also admitted to the second nozzle. Finally 
under periods of near full load conditions, steam is 
admitted to all three nozzles. All of the steam admitted 
to the first nozzle also passes through the second and 
third nozzles. The nozzles convert the potential energy 
of the steam to kinetic energy by accelerating the vapor 
molecules. Kinetic energy of the vapor molecules is 
transferred to the rotating energy of the turbine by 
Steiming the turbines blades. The first two stages 
(after the first nozzle) receive the steam at its high- 
est temperature and pressure and convert the energy with- 
out allowing the steam to expand further. These stages 


are called impulse stages. They reverse the flow of vapor 
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molecules without changing its velocity relative to 

the rotating blades. The remaining seven stages (three 
after the second nozzle and four after the third nozzle) 
are reaction stages. These blades also reverse the flow 
of steam to recieve its energy, but additionally they 
act as partial nozzles accelerating the vapor molecules 
as they pass. This 1S accomplished by the exhaust Side 
of the blades having a larger opening than the inlet 
Side. The resulting drop in presSure accelerates the 
steam. 

The majority of the steam passes through the first 
nozzle and down the entire length of the turbine. There 
is a delay from the time the steam is admitted to the 
steam chest around the tirst nozzle until at actuallv 
passes through the nozzle. As the steam reaches the 
turbine blading it imparts more of itsS energy at each 
Stage. The turbine power iS measured to be linearly 
proportional tothe steam flow. (See graph 1.) The turbine 
time constant is associated with the delay and transient 
time of the steam passing through the turbine. The 


turbine time constant, ue controls the rate of change 


h’ 
of the turbine power, TP, caused by the throttle position, 
VP, [3]. 

TP = (VP - TP)/Top 222 


By integrating this equation instantaneous turbine power 


may be found. 
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Pee ida MECHANICAL SYSTEM 

From the current turbine power a value for the 
speed of the turbine may be calculated. Since power 
must be conserved the difference of the turbine power 
produced minus the electrical power and the lost power 
must appear as a change in angular momentum. The only 
method of changing the angular momentum of the turbogenera- 
tor is by accelerating it. The acceleration of the 
turbogenerator, 0, is equal to the difference of the 
power produced and the power required divided by the 
inertia constant, 2H. This acceleration is also the 
acceleration of the power angle, 8. Since Torque times 
speed is equal to the power: 
ee = ( TP — f te + TK. 1* «Zon Zee 

The inertia constant, H, is computed from the 
moment of inertia of the turbogenerator and its rated 
speed. By knowing the dimensions and weight of the 
rotating parts the moment of inertia can be computed. 
However, the difference in the speed of rotation between 
the turbine and generator must be taken into consideration. 
The larger radius of the generator and gear dominate the 
inertia constant despite the larger speed of the turbine 
Sane plm1oOn . 


The term in equation 2.3.1 is an approxima- 


Tioss 


tion of all the mechanical losses. The largest part is 


Gue to the speed reducing gear. Friction losses in the 





bearings are also included, as well as windage losses 
in the generator. The sum of the losses was assumed to 
be equal to the total power required to maintain rated 
speed with no field excitation. This value is Th oss 
weed, Or 1000 lbs of steam per hour, [2] .- 
Poe tHE SPEED CONTROLLER 

The throttle position is controlled by the governor. 
The Watt governor measures the turbine shaft speed by 
means of a flyball arrangement. Attached to the shaft 
of the governor are two opposing balls on a pivoted 
linkage. The balls may move easily in a radial direction 
and are forced to rotate with the shaft. The two balls 
are connected by a spring. The spring exerts the 
rorce regquirea to make tne palls travel in a circular 
path about the shaft. As the shaft's speed increases 
more force 1s required to continue to accelerate the 
balls to their circular path, and the spring stretches 
increasing the force it exerts, and enlarging the radius 
of the circular path of the balls. Another linkage 
connected to the balls detects the increase in the radius 
of the circular path of the balls and uses this movement 
to displace a valve. The movement of the valve causes 
a pressure difference between the hydraulic pressure in 
the power piston tending to open the throttle, and the 
SoraingeLLying to close the enroeeie: The power Sister 


moves to equalize tne forces and in so moving returns the 
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valve to its neutral position. The rate of change of 
the throttle depends on an error signal from the 
governor, l-w, amplified by the hydraulic system and 
combined with the old throttle position, DVP, minus the 
current throttle position, VP, divided by the time 
constant of the hydraulic servo system. The delayed 
throttle position may be found by integrating the 
throttle position minus the delayed throttle position 
divided by the time constant, [3]. 

VP = ( DVP - VP + [ l-w 1*35) /., , 2.4.1 

DVP = ( VP - DVP )/305 2 wAene 

The governor maintains speed within 4% of rated 
value and returns to a steady state speed of within 23 
of rated value in less than two seconds under normal 
Boag mluctuations. 
2.9 THE GENERATOR 

The Allis Chalmers generator is a 1000 KW, 0.8 
Power Factor, 450 Volt 3 Phase, 60 cycle totally 
enclosed synchronous A.C. generator, [2] and [7]. It is 
Seaeltent pole type machine operating at 1200 RPM. The 
rotor receives direct current through four brushes on 
two rings from a static excitor. The rotor has complete 
damper windings in the pole faces. The stator has 90 
coils, connected in groups of five coils. There are six 
of these groups per phase with the phases connected in 


a Y circuit (see figure 4). The air gap is approximately 





1K: 


0.185 inches. Both the rotor and stator are punched 
from electrical sheet steel with each sheet insulated. 
The generator 1S air cooled. 

Three basic asSumptions are made in modeling the 
electric and magnetic circuits. The first assumption 
is that each phase 1s identical and that each coil in 
each phase 1s identical so that the model will be 
described as a two pole machine (see figure 5). 
Even though the amortissieurs are continuous about the 
rotor, mathematically they will be viewed as being 
symmetric about the direct and quadrature axes so that 
no mutual components exist between the direct and 
quadrature axis amortissieur. The final assumption, in 
steady state, all Lf1ux distributions vary aps siiusuirus 
momespace, (4) and [5]. 

According to Faraday's Law, the voltage induced in 
a circuit iS equal to the rate of change of flux 
linkage. Assuming the resistance of the three phases are 
equal, the phase voltages, e ep © are equal to the 


a eld 


meee Of change of flux, d*/dt, minus the phase current, 
i 


ae 2 times the phase resistance, r. 


a cy 
oon = tH ay Aene - gh 
en = d¥p/dt - rip } 2c 


The flux linkages in the phases result from the currents 


Bilevimag in all the coils of the generator. The self and 
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1b, 


mutual inductances of these coils may vary with the 

rotor angle, 8. Since the three phases are assumed to be 
identical except for being 120 degrees apart, only the 

a phase will be described. By adding 120 degrees to 

the rotor angle at every occurence the inductances of 

the b phase will be found. Likewise by subtracting 120 
degrees from the rotor angle the inductances of the c 
phase will be found. The self inductances of the a phase 
is expressed by an even Fourier series of the rotor 


angle. The first term, a constant L results from the 


aaQ’ 
actual inductances of the a phase coil by itself. The 
Poeeecerm, Li _5cos20, is proportional to the cosine of 
twice the rotor angle. This term is solely due to the 
Salience of the rotor alternately decreasing and increas- 
ing the reluctance to the stator coils. Since the sum 

of the frequency components, other than the component at 
60 Hertz, of the terminal voltage and current 1S measured 
to be less than 2% of rated value, only these first two 
terms of the Fourier series will be considered. If the 
frequency components of the generator output were to be 
Studied then the Fourier series would have to be extended 
to encompass more terms [6]. The self inductance of the 
a phase is -L,.9-1,,7¢0829. It is negative so that 
current out of the generator will represent positive 
power. (See figure 5 for all sign convention.) The 


mutual inductances between the phases is also expressed 





as a Fourier series. The constant term, Labo: is solely 
due to the phase coils alone. The same second term 
appears in the Sap iee. However the reluctance is now at 

a minimum when the pole face is half way between the two 
phases, a to b Thee + Loq2Cos28t+60; a to ¢ Lang + Laa2COS 
28-60. The mutual inductances between the rotor and phase 
a Gepend completely on the rotor angle. Since the field 
coil and mathematical direct axis amortissieur act along 
the same axis they appear mathematically similar. The 
Mutual inductance from the field, L,fcosé, and from the 
direct axis amortissieur, Lanqgcos®, are at a maximum 

when the north pole, positive flux, of the rotor is inline 
with the a phase. The quadrature axis amortissieur 
wenaves Similarly, eacepe that Lt 1S Mathemacicaily 720 
degrees ahead of the direct axis. The mutual inductance 


of the quadrature axis amortissieur, L sin@, 1S a 


akg 
maximum 90 degrees before the a phase is inline with the 
rotor. Again the inductances of phases b and c are 120 
degrees ahead and behind the a phase. Since the stator 
is totally symmetric (disregarding slots) the self and 
Muteal inductances of the rotor coils are constants, Lr 
mommene tield coil; i,q for the direct axis amortissieur; 
and Lig for the quadrature axis amortissieur. There is 
mutual inductance between the field and direct axis, 


Ee Pue Since the quadrature axis is 90 degrees apart 


there is no flux linkage between them. From this 





Zi 


maeOrmatLoOn an inductance matrix, L, may be formed. 
The flux linkage, ¥, is then equal to the inductance 
times the current, or in matrix form: 
y=eLFi 2. ee 

Pececigure 6 for Lh. 

Biceehase voltages, €4 €, €,, Nave’ been discussed 
once in equation 2.5.1, however they may also be found 
in relation to the generator's load. That is the 


wemedadesia equal to the current generated, 1, 1, lo, 


times the load, Rigad © Sioad- 


= * 4 
Sa See ta : Ly oad 


se de ip ya CL Meh AS 8 


s di/dt 
Cp = Rioaa * 
Qc = Rioaa * ic + Loaa * dig /dt 

Tuc field Voltaye a5 Ucsaveu FLOm Cue cacactus’ os vulpur 
and will be discussed later. The amortissieur coils 
are shorted. 


The generator may be modelled by the following 


Matrix equations [4], [5], and [8]: 


aydt = e — r™i 2 oe 
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diyvdt = 


es 


ZS 


(L + Lygag)) *C @ + E*L - SR/dt *i) 2.5.8 


This three phase model may be simplified by the 


use of a coordinate transformation. The new equations, 


known as Parks' 


EO equations 2.5.1 to 


Ze Deo 


equations, will form a model identical 


The coordinate transforma- 


tion will remove the dependance on rotor angle from the 


BAauetanece Matrix Lb. This model will be called the dq0 


model. The transformation is as follows: 
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The current, voltage, and flux linkage matrices for the 


dqO model are found by: 


7 =e eee 
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s el 
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The inductance matrix for the dq0Q model is found by: 
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All other quanities are 3/2 their three phase values. 
This model is advantaceous because it is time invariant 
however it is still nonlinear because of the cross 


eebms OL frequency and flux. 


Ve. = 5 9 oe * ig + uw * Yq 
mee = &g - = * eam 2.5.15 
dy = < 

O/ at = € iets 16 


The torque of the generator may be derived from 
the rate of change of air gap power [41 or for the 
three phase model 
Tere = “/37gl¥alicnip) t¥p (ig-ic) +¥o(ip-ig)] 2.5.16 
For the dgq0O model 
Tele . “a 7 GMa 2 - ule! 
Zeooee tHE VOLTAGE CONTROLLER 
The excitor and voltage regulator are combined into 
one unit in‘this system. In all but the initial starting 
of the turbogenerator, the excitation is derived from 
the rectified terminal voltage and current. The initial 
field flashing voltage is obtained from a permanent 
magnet alternator attached to the generator shaft. At 
times other than starting it is used as a tachometer to 
the main switchboard. It is capable of supplying enough 
excitation to produce one half the rated terminal 


voltage at no load and is disconnected electrically 


from the system at all other times. 


[Oo 
an 
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The normal excitation is derived from a set of 
three magnetic amplifiers, one for each phase. The 
magnetic amplifiers have three inputs, a saturation 
input, a Signal input, and a control input. A current 
transformer AGTOCS ONG Phasecwer the Generator provides 
two of these. The current transformer has windings 
both parallel to and in series with the terminal. The 
potential windings Saeuide power when the generator is 
in a no load condition, whereas, the current coil provide 
power if a short circuit develops. Under normal conditions 
the current transformer produces an output proportional 
to the vector sum of the terminal voltage and the terminal 
current. The transformer's output is rectified to provide 
emer saturation imew. Mee tie magnet amel¥i see. ive 
output is also the signal input to the magnetic amplifier. 
The control signal is derived by an error detecting 
circuit to be described later. The output of the magnetic 
amplifier goes to a three phase rectifier with the other 
two magnetic amplifiers. This rectified output is the 
source of excitation voltage. There is a free wheeling 
diode across the field terminals to provide a path for 
the field current during short periods of low voltage 
from the magnetic amplifiers. 


The actual control of the magnetic amplifiers is 


= 
—_ ae 
= 
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Gerived from a pair of potential transformers which 
form an open bridge in the generator's terminal lines. 
The output of the transformers is input to a voltage 
sensing bridge of two zener diodes and two resistors. 
The voltage error signal 1s amplified by a magnetic 
amplifier. A damping signal, derived from a RC 
aifferentating network in the excitor terminals, is 
subtracted from the error signal. The resulting 
Signal is smoothed and used to alter the gain of the 
three main magnetic amplifiers. 

Two separate actions occur when the load of the 
generator is altered. First the input signal of the 
main field magnetic amplifier is modified. AW olalis 
as the vector sum ol terititnaL VOLtaye ana Current. 

As the load is increased the terminal current increases, 

increasing the vector sum. If the terminal voltage 

were to remain at rated value this change alone would 

maintain the necessary excitation. But the terminal 

voltage does not stay constant. As it changes the 

error detector circuit alters the gain of the magnetic 

amplifier to try to return the terminal voltage to 

rated value. These two signals working together improve 

the response of the terminal voltage beyond the 

response either signal alone is capable of producing. 
The model for this static excitor and voltage 


regulator is based on a IEEE Standard type three system 
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meomero). The Output of the current transformer, TER, 
1s the sum of the magnitude of terminal voltage, VT, 
plus the magnitude of the terminal current, TI, times 
the power angle. Since the values are rectified the 
actual phase angles are not important. The error 
detector circuit produces a Signal, ERR, which is the 
rated per unit terminal voltage minus the normalized 
terminal voltage. The stabilization signal, STAB, is 
Bemnd trom the derivative of the field voltage, ef, 
integrated with the old stabilization signal. The 
change of the magnetic amplifier gain is integrated 
from the old gain, SIG, and the error signal and the 
Stabilization signal. Thus 


aes Ex. Vile =~ Pe 


an = 9 lb =aV i 


e K i. 
SOE =! STAB GGe7 * SF BS ail 
see =) 7 * SIG + Ka/, * (ERR-STAB) 
a a 
e, = (1+ SIG) * TER 


(See figure 7 .) 
Ka and Ta are parameters of the main field magnetic 
amplifiers. Kf and Tf£ are parameters of the error 


Grreulte Magnetic amplifier. 
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ae) eo weeeETON OF THE MODEL 

The models were assembled by parts to form a 
complete model of the turbogenerator. The model was 
Simulated on a digital computer by a fourth order Runga 
Kutta type algorithm. It was not a true Runga Kutta 
scheme, Since the solution was found in two steps. First 
the generator model was solved for one step under constant 
excitation and rotation. The voltage regulator, the 
speed controller, and the mechanical system were stepped 
forward one step under constant generator output. Each 
step provided the input for the next step. The method 
provided satisfactory simulation because of the large 
differences in the time constants of the two parts of the 
Syscem. In this manner wne two genéiator nedeis Cuuid 
be easily interchanged as well as the new controllers. 

The three phase generator model required the 
smallest step size and the largest computing time. 

The smaller step size resulted from the currents and 
voltages in this model remaining Sinusiodal. The longer 
computing time resulted from the inversion of the 
induction matrix at each step. 

The advantages of the three phase model were that 
the outputs are actual machine quantities. The model 
would also accept any conceivable load conditions. The 
model is capable of keing more accurate in its simula- 


tion since exact coil inductances may be included. This 
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allows the model to accept possible differences in the 
three phases. 

The dqOQ model can be simulated in much the same 
manner as the three phase. However, the dqO is time 
invariant, thus the inductance matrix need only be 
inverted once. By proper matrix manipulation the 
required computations for this model may be reduced. 
Since the transformation removes the frequency 
component associated with the rotor angle ( the 
highest frequency oscillation and the only undamped 
oscillation) a varible step size may be used to increase 
the step size as the oscillations disappear. All of 
this results in a model that the computer may execute 
Mees Laster than tee three phase model. 

The dqO model is restricted in simulating unbalanced 
loads. To accomplish this the output must be converted 
to phase quantities, the load values computed and then 
transformed back to dg0O values. This still requires less 
time to accomplish than the three phase model requires 
to invert the inductance matrix. The dqO model may not 
include any differerices between the three phases. 

In future studies saturation should be included in 
both models. In the three phase model the inductances 
of each coil would be varied by the current flowing in 
them to represent the saturation. This would not 


significantly affect the model. However, in the dq0O 
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model if the inductances had to be altered for satura- 
tion the matrix would have to be reinverted at each 
change in saturation. This could significantly reduce 
the speed of execution of the dq0O model. 
2.8 TESTING THE MODEL 
Both models were tested under the following 

Gonditions: 

Nemoad to fully load 

NoOmtOoad "CO. SROGE C1 LeulLt 

Full Load to partial load 
The models were found to compare satisfactorily in all 
respects with each other as well as with the provided 
acceptance tests for the turbogenerator. The actual 


Simulation run results may pe found in section 4. 





TIf THE DEVELOPMENT OF THE DIGITAL CONTROLLERS 


In the introduction several proprieties of the ship- 
board environment were introduced which would enchance 
a simple control scheme. The elimination of the boiler 
from the model has been discussed. This also allows the 
Speed controller to operate without determining the 
steam power available. It appears to be a constant 
infinite supply. The voltage regulator will not have to 
be able to dampen oscillations caused by a multi-machine 
environment since there appears to only be one machine 
and one load. Finally both controllers have leeway in 
the control. The output does not have to be exact at 
all times. 

ieee ten yn ene eaameempe +o bul ld eerese 
controllers. Only the simple control scheme will be 
proposed. Developing the hardware required to actually 
build the controller will be left to a future study. 
3.1 THE DIGITAL VOLTAGE REGULATOR 

The first and most obvious method of finding a 
digital control scheme is to implement the present 
controller in digital form. The present controller 
consists of an error detecting circuit. The main field 
magnetic amplifier provides the power source for excita- 
fone Tie proposed digital controller will only control 
ne woutput of the magnetic amplifier. The current 


controller model is: 


= w 





STAB = -1/p, * STAB + Ke/pe * &¢ 25633 
STG = -l/q, * SIG + Ka/p, * (ERR-STAB) 2.6.4 


In digital, z transform, form with a sampling period of 


T the following transfer equations would result: 


STAB(z) _ Kf/Te Z, 


SIG(z) Ka/Ta z 
(PEp=sirAB) (2) > 5 2 e-2/T, Sale ge 


The difference equations provide a better feel for what 
this proposed AVR would require: 

STAB(K+1)T] = Ke/re*e-[(k+1)T]t+e7'/T£ *STAB[kT] 

SIG({k+1)T] = Ka/p_*(ERR-STAB)[(k+1)T]+e7/Ta SIG{kT] 
To implement this proposed digital information on AVR 
the terminal voltage, or error signal, and on the field 
voltage, or its derivitive, would have to be provided 
digitally. The former signal in each case would require 
processing to arrive at the needed signal. The voltage 
regulator could be greatly improved by the use of some 
functions easily performed digitally but which in the 
analogue world would be extremely difficult. A very 
Simple generator model would be, [15] 

ce =) EX Beis Ss 

Where VI is terminal voltage, I is terminal current, X, is 


synchronous reactance and k is a constant. This model 
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- Shows that the required excitation current is related 
to the real current. Thus if the real current could be 
calculated the exact excitation current would be known. 
The real current depends on the exact load and the rotor 
poe Lt 1S difficult to compute. 

An eaSier method of approximating the required 
excitation comes from equation 3.1.3. By measuring the 
instantaneous terminal voltage and field current an 
approximation of the field current required to return to 
the rated terminal voltage conditions would be eo 
This assumes that everything else remains constant. To 
see this the model must be examined. By normalizing the 
rated terminal voltage to one then forming a proportionality: 

VT a iceewe age Xe 


ay = “Kedh — 4 T'#X Serle 
& Ss 


(The prime quantities are the desired values.) 


1 
But under constant rotor angle: T' = [se 


a 
Thus the desired relationship: if = ier 


Although this is a good approximation, it is not 
exact. The major fault is that neither the rotor angle, 
nor the load, remains constant. The real current 1s 
therefore not always SrOnomuon ey to the terminal 
moltage. However, since the field current cannot be 


changed instantaneously, several new approximations can 





be made updating the required field current before the 
generator can completely readjust to the new load. 

To insure this procedure works correctly it must be 
Shown that it quickly produces the correct estimate of 
the required field current, and that it will maintain 
this correct estimate. The estimate can be viewed as 
being proportional to the instantaneous power demanded 
of the generator. Under this instantaneous power 
requirement the excitation current must be increased to 
increase the terminal voltage. In some possible load 
changes (such as application of a large resistive load 
from noload conditions) the instantaneous power re- 
quirement is not even in the same direction from the 
previous power required as the steady state power 
reguarement of the mew ivau. (Under appiscacion of russe 
resistive load the intial terminal voltage rises 
indicating less power is required although the steady 
State power requirement is larger than the present 
requirement.) The estimate of the required excitation 
must be updated regularly. At the next estimate the 
instantaneous power requirement and rotor angle will be 
closer to the steady state values, and will produce a 
better estimate to the required steady state excitation. 
Requiring a new estimate to be calculated before the 
field current can significantly change will insure that the 


power will not be changed in the wrong direction. 
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The procedure might be better explained graphically 
in figure 8. Here the generator is at steady state under 
the original power requirement. The load is then altered. 
famee the field current cannot change instantaneously the 
terminal voltage drops toward the new load line. The AVR 
assumes the load has not changed and raises the excitation 
voltage proportional to the terminal voltage, to return 
the terminal voltage under the original load to its rated 
value. At the next approximation the field current is 
increased and the terminal voltage is still below rated 
value. Notice that the approximations are at a high enough 
rate to prevent the field current from Significantly chang- 
ing between approximations. Even though the field current has 
been increased the terminal voltage 1S not at its rated value. 
The AVR estimates the new power requirement and computes a 
new estimate of the required field current. The process 
converges rapidly to the correct terminal voltage, and the 
AVR no longer attempts to change the field current. The AVR 
thus always produces a better estimate of the required field 
Current and steadies out at rated value. 

Although this proposed AVR works, it is slow to 
respond because of the large time constant of the field 
Gurenic. To ee the rate of change of field current 
the excitation voltage applied to the field must be larger 
than required. An easy method to overcompensate the field 


voltage is to use *£/ (yer)? as the voltage signal. 





/new excitation 
requirement 






ee ee 


- second est. new 
excitation req. 





tem R Ae O (0) 
}~ 


£2 
“£5 first est. new 
1¢ i ave tatiryy req. 
| . 
pall original 
ie excitation req. 


VT VT, 1 


Terminal Voltage 


( The * represent new approximations are being made by the AVR. ) 


Figure 8 
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Then the voltage will be exact when at unity and will 
be overcompensating at all other times. The field 
current will therefore respond to change faster than 
if just the required voltage were applied. To improve 
even further, if the terminal voltage varies more 
than a preset amount from rated value the AVR would 
push the excitor to its limits. Together these changes 
make the proposed digital AVR respond much faster, but 
the cost is that the terminal voltage tends to oscillate. 
To prevent this the AVR must calculate a new estimate 
of the required field current before the field current 
can significantly change under the maximum rate of 
change of the excitor. The samples must be often 
enough to prevent the terminal voltage from changing 
more than this preset amount before another sample is 
taken. The excitor voltage can then be set to a better 
estimate of the required excitation before the field 
Current can pass through its required level. A 
sampling period of once every three radians, twice a 
cycle, appears under most conditions to be an adaquate 
aise . 

To implement the proposed digital automatic voltage 
mequbator, the following tests of the terminal voltage 
would be made: 


maximum 


VT 


1A 


0.97 then Cr 
Vi wa ae then er = minimum 


1 


VT = 0.0 then Cr 
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If these tests fail, i.e. VT is in between .97 and 1.03 
then the excitation voltage, e¢, would be equal to the 
mcler current, le, divided by tne terminal voltage 
squared: 
e¢l(k+1)T] = i, {kT] / (VT[kT])? 
The test for zero terminal voltage is to insure that 
meer Short circuit conditions there is a field current 
maintained so that the short circuit current will last 
long enough to trip the circuit breakers about the short, 
[10]. A comparison between this proposed digital AVR 
and the present controller can be found in section four. 
3.2 THE PROPOSED DIGITAL SPEED CONTROLLER 
A digital governor could also be found by implement- 

ing the present controller in a digital scheme. This 
digital governor would require a digital tachometer or 
frequency counter and a digital position indicator for 
fiemeirottie. The linearized model of the current 
governor is: 

VP = (DVP - VP + (1- »)*K,)/T, Deedee 

DVP = (VP - DvP) /T, Diane 
These equations would convert to the following z transform 


transfer equation: 





VP (z) 
Gia) (Zz) ~ _ TT T) es 
22 - [1-¢ pt} e+e HIN 
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In difference equation form 


VeIGk+1)T)] = (Ite 13) vp [kT] - a 8yp(k-1) 7] = 3 ‘h/T,, (1-w) (+1) 7] 


a ~T8 
Ti (Thtl1) > Th(TrtTD (1-w){ kT] Sa2c2 


A much simpler digital governor can be found. 

The present governor depends on a speed change of 
fixed magnitude before the throttle reacts to a change 
in generator load. Since the inerita constant of the 
BuEpegeneraror 15 SO large this governor reacts 
sluggishly to load changes. The result is that the 
turbogenerator is allowed to oscillate slowly about the 
desired frequency. A digital governor could easily 
improve this action by altering the turbine power 
peemeracithe speed hats = ohteee to chit@e. Tail  peereecd 
digital governor would compute the power output of the 
generator and use this information to alter the turbines 
power. This would greatly reduce the speed fluctuations. 
The easiest method of arriving at generator power is the 
instantaneous generator voltage times current. £4) 
complete the governor an actual speed signal must be 
used to maintain the exact speed. Thus an error signal 
can be added to the generator power to return the system 
to rated speed. However, when an error signal is added, 
the speed will again fluctuate. To prevent this, 
information about the change of speed of the system can 


be subtracted from the error Signal. In this manner the 





speed can be returned to rated value and not oscillate. 


The proposed governor would use generator power, E Salar 


and system speed, »w, to control throttle position, VP: 


VP[(k+1) T] 


{Pozo (kT] + -1 + (?/qr-l) wI[kT]-°/qul(k-1)T] 3/top, 


, VikT] ~2ese) 


Pore [kT] 


A comparison between the proposed digital governor and 


the present governor can be found in Section four. 
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iV. THE RESUS TOS Ine oul ee ONS 


In this section the models and controllers will be 
compared and evaluated. First the models will be 
compared to the actual acceptance test results, and to 
Seaciwotner. Then the proposed digital controller will 
be compared to the present controller. Finally the 
remaining test results will be presented. Most of the 
information will be presented graphically as the 
responses of the models to various imposed loads. 
fee lLsolto ON THE MODELS 

The generator acceptance tests provide information 
on the actual generator response to three phase short 
empeuit.anad to nOload to full load, full volt-amperes 
lige. LRaS suit ivau was Gn szuduction motor witn its 
rotor locked. It provided a per unit impedence of 
.5887 + 31.691. The generators response May be seen in 
figures 4.1.1 and 4.1.2. The tests were on the generator 
alone and no turbogenerator test results were available. 
The frequency responses therefore are based on the design 
specifications. In figures 4.1.3 and 4.1.4 the dq0 
generator model was simulated under the same load 
conditions as the acceptance tests. Under the no load to 
full load the dq0 model produced smoother results but 
they were similar to the test response. Table 4.1.1 


provides a summary of the responses. 
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Tnese tests were repeated and the results are illustrated in 
figures 4.1.5 and 4.1.6 with the three phase model providing 
the responses. Again the resultS are similar. Because of the 
amount of time required to simulate this model, the responses 
are for shorter periods of time. The accuracy of the model 
1S apparent from these results. 

Comparison of the two models to determine which is 
more accurate is difficult. Both models perform 

adequately under the information available for these tests. 
The determining factor in choosing which model to use 
should be based on the advantages of the specific model 
to Simulate the conditions which will be imposed. The 
dgQ model has advantages in speed of execution and 
memory required. The dq0 model requires on the average 
0.066 CPU seconds per cycle to Simulate average load 
conditions whereas, the three pnase model required 1.6 CPU 
seconds to simulate the same load. Also the dq0O model 
required about two thirds the memory that the three phase 
model required. Under other load conditions the difference 
in execution time can be reduced but the three phase model 
always requires more time. 
eee tests ON THE PROPOSED CONTROLLERS 

Now the simulation of the proposed digital controllers 
will be compared to the model of the present controllers. 
Because of its higher speed of execution the dq0O model was 


Msedethroughout this section. First the Erequency response 
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generator -0185 42 = fe 
dqO model 2023 . 333 .70 
three phase model 7020 :333 XXX 


Table 4.2.1 
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St the turbogenerator to the no load to full load, locked 
rotor induction motor, must be established. The design 
specifications is for the turbogenerator to remain within 
6% of rated value, 4% for normal load fluctuation, ane 
return to within two percent of rated speed within two 
seconds. Figure 4.2.1 is the simulated response of the 
system frequency under this load. 

In figure 4.2.2 and 4.2.3 the dqQ generator model 
is controlled by the proposed digital controller. 

Again the voltage dip is .89, however the rise time and 
settling time are greatly reduced. The frequency 
response is even better. Here the variation is one third 
as great as the actual controller. The proposed digital 
controller also returns the frequency to exactiv rated 
value, rather than to within two percent of rated value. 
The proposed governor requires about three seconds to 
return to rated speed. 

To further test the proposed digital controller an 
even larger load fluctuation was simulated. The locked 
rotor induction motor provided full volt-amperes for the 
generator but the actual power required was only about 
forty percent of the available power. In this test both 
the volt-amperes and the power required were a maximum. 
The simulated load was .8 + j .6 per unit. The simulated 
responses under the present controller are shown in figures 


4.2.4 and 4.2.5. The same load conditions were repeated 
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with the proposed digital controller in figure 4.2.6 
and 4.2.7. The proposed digital controller was faster 
in its response to load changes. 

Finally to insure that the proposed digital 
controller worked completely the controller was tested 
under a decrease in load: from a full power load the 
locked rotor induction motor load. These test results, 
Figures 4.2.8 and 4.2.9, were important because the 
transients induced in the field coil make the proposed 
AVR overshoot the correct value. The proposed digital 
controller still performed well on this test. fThis 
Simulation was carried far enough to be able to see the 
governor return the system to rated speed. 

Threanehant thie and ald other tests the cimnle 
proposed digital controller out performed the actual 
machine controllers. 

4.3 TESTS ON SAMPLING RATE 

The rapidity at which the proposed digital 
controller must reestimate the correct throttle position 
and excitation voltage greatly affect the cost of the 
proposed controller. For the digital automatic voltage 
regulator the full resistive load is the worst case. 
Under this load the generator responds faster than any 
other condition. By simulation under full resistive 
load if the digital AVR reestimates a new required 


excitation voltage every 0.4 cycles, a new estimate of 
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the required excitation voltage can be calculated before 
the terminal voltage can overshoot its rated value. 
Thus, the field current does not overshoot the predicted 
value even when the terminal voltage is below the value 
forcing the excitor into its limits. 

Under other load conditions, especially high 
inductive loads the period between estimates can be 
considerably slower. Under the locked rotor induction 
motor load, the proposed digital AVR would only have to 
reestimate the required excitation voltage once every 2.3 
cycles to prevent the terminal voltage from oscillating. 
For the load .8 + 3.6 pu. the AVR must reestimate the 


required excitation voltage at least once a cycle to 
nrevent the terminal wonltace from oscillating. To 
insure that the AVR does noe cause oscillation of the 
terminal voltage a sampling rate of 120 times per second 
1s recommended. This rate would adequately handle any 
load condition. 

The proposed digital governor does not have such a 
critical sampling rate. The sampling rate of the governor 
would depend more on the frequency of the load changes. 
The slower the sampling rate the larger would be the 
frequency fluctuations. A rate of about thirty times a 
second or once every two cycles is recommended to allow 


this proposed digital governor to out perform the — 


present governor. 
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4.4 POSSIBLE FAULTS IN THE PROPOSED CONTROLLERS 

Another problem that might arise with the proposed 
digital AVR could result from not correctly knowing the 
field resistance. This would cause the AVR to incorrectly 
estimate the required excitation voltage. This would 
result in either an oscillation if the field resistance 
were lower than expected or in the AVR's response time 
decreasing significantly if the resistance were larger 
than expected. The error in the value of the field 
resistance would have to be considerable to noticeably 
affect the operation of the proposed AVR. If the 
terminal voltage did tend to oscillate slightly and 
this was determined to be the cause, the AVR could be 
altered to use information about the rate of change of 
field current to improve its estimate of required 


excitation voltage and eliminate the oscillation. 





V CONCLUSIONS 


The objectiveof this thesis was to compare the three 
phase generator model to the dqO generator model and to 
find a simple digital control scheme which would meet 
the design specifications of the shipboard turbogenerator. 

The two computer models, the three phase and the 
dq0, were found to agree completely in all respects. 

The three phase model would require more computer storage 
and execution time than the dqO0 model. However, the 
three phase model would simulate more complicated load 
and saturation conditions than would the dg0O model. 

This thesis showed that if complicated load and 
Saturation conditions were to be simulated the three phase 
mOGcl would be boviericiul. if baianced load aiu satura- 
tion conditions were to be simulated the dqQ model would 
result in considerable savings in time. Botn models 
were satisfactory in simulating the different controllers. 

A simple digital controller has been proposed which 
by using information easily obtained digitally, but 
difficult to obtain by analogue means, has significantly 
improved the performance of the turbogenerator model to 
load changes. 

This investigation should be continued to further 
test the affects of saturation and a multi-machine 
environment. After these tests were satisfactorily 
passed the supreme test wouid be to build and install 


Swemrvarndagital controller. 
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APPENDIX 


weeecalcvlation Of Inertia Constance 

To calculate the inertia constant, H, the moment 
of inertia of the system was required, Figure 9 is a 
basic drawing of the rotating parts showing size and 
weight of each. By assuming constant density through- 
Out each piece and neglecting light projections, such as 
the coupling and actual turbine blades, the moment of 


inertia was computed by integrating 


L Jaa (AS 
WR2 = | | o r? drdéedl 
0 0 0 


Or Since all parts are cylinderical WR? = mr /, 


under constant density. 


The moment of inertie of the matanr narte ic: 


turbine 89.08 ft ps at 10009 RPM 
em | ~89 fits 1s at 10009 RPM 
gear 81843.8 ft ibs ae 2 U0 REM 
generator 8680.8 ft ibs aiteyee2 Cer 


The inertia constant is the moment of inertia converted 
to kinetic energy terms and put in the per unit evecare 
Thus converting to megajoule-seconds per radian and 
then dividing by the apparent power, volt-amperes, 


the inertia constant is, |{11]: 


H= 24.29 so aeeeee = (actually dimensionless) 
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A.2 Calculation of Generator Constants 

wie shert circuit erests lie prevuded the 
information to calculate the machine constants. The 
Short circuit provided values for direct axis reactance, 
Kar damect: axis transtentereaceance - Xa direct axis 
Subtransient reactance, Xe and external reactance, Xo 
Also given were the phase and rotor resistances and the 
voltage dip for an initial load. From this the fact 
that x, 7 x3 would indicate that the generator was a 
Salient pole machine which was previously known. The 
subtansient quadrature axis reactance was calculated 
from Xo and Xa by Xo = Xa XG . The external reactance 
is an approximate time average of the subtransient 
reactances. The quadrature axis reactance was more 
difficult. The only information available was the 
voltage dip, however this depends mainly on the direct 
axils transient reactance. This was used as an 
approximation for the value of the quadrature axis 
reactance, Xq° The zero axis reactance, Xs which was 
never used since only balanced loads were simulated was 
derived from Xo Xar and Xe 

With the machine time constants the remaining 
constants followed easily. The direct axis to field 
axis mutual inductance was assumed to be only slightly 
less than the direct axiorneaseanes due to the leakage 


of lines of flux. The same leakage was assumed for the 


quadrature axis. Thus 
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xX = 


The field time constant, T! , by 





d 
Le xq 
ie = 139 
2 Re Xq 
The field transient time constant, Ts Dy. 
2 
aa L r $8 
Ly £kd/y, sl 
ae _ —— 
4 Bo x! 140 


The quadrature axis time constant had to be completely 
approximated by a table in [4], chosing a value for a 


machine which was Similar in the known constants. 


Z 
L 
Lkg - GR/1, 


T" = 138 
e | Reg 


The three phase constants are two thirds of these values 


except for the stator values. These phase values are 


determined by 
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= 3 
Xa = Laao + Lapo + “/o2baap ze 
- sg 
a Load i Labo /yusar 26 
6 ~ eae - 2 eee 27 


These equations are directly from [4], and [5], the 


equation numbers are the numbers assigned in [4]. 


Machine Constants 


Measured Values (again all values are in per unit ) 


Ke = ie | = 1078 radians 
@ = " = : 
Xa 0.167 v3 7.78 radians 
x4 = 0.135 Voltage Dip = .89 
= 4 
Xo Ca23 


Computed Values 


dq0O Model Three Phase Model 
L = sD Lad en) e916 

i = 1.217 Lb = 0.4225 

Ly = 0.091 Loa2 OF 0943 

Lae = 1.4 I =) 5 BIRyE 

Laxq = 1-4 Lakq = -933 

Lokq = 1-14 Lakq = -76 

L¢ = 1.47 Lf = .98 

Lkq =1.47 Lkq = .98 

Lkq = 1.217 Ups =) eile: 


aos Lfkd oso / 


Lfkd 





Rea 


I! 


UPOuIs 725 
O20 Fyne 
0.143 


OF e 8)aes 
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0.000915 
0.01142 
O09 15 
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